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ABSTRACT 

We present the discovery that UScoCTIO 5, a known spectroscopic binary in the Upper Scorpius 
star-forming region (P = 34 days, Mjot sin(i) = O.64M0), is an eclipsing system with both primary 
and secondary eclipses apparent in K2 light curves obtained during Campaign 2. We have simul¬ 
taneously lit the eclipse profiles from the K2 light curves and the existing RV data to demonstrate 
that UScoCTIO 5 consists of a pair of nearly identical M4.5 stars with Ma = 0.329 ± O.OO2M0, 

Ra = 0.834 ± O.OO6i?0, Mb = 0.317 ± O.OO2M0, and Rb = 0.810 ± O.OO6i?0. The radii are broadly 
consistent with pre-main sequence ages predicted by stellar evolutionary models, but none agree to 
within the uncertainties. All models predict systematically incorrect masses at the 25-50% level for the 
HR diagram position of these mid-M dwarfs, suggesting significant modifications to mass-dependent 
outcomes of star and planet formation. The form of the discrepancy for most model sets is not that 
they predict luminosities that are too low, but rather that they predict temperatures that are too high, 
suggesting that the models do not fully encompass the physics of energy transport (via convection 
and/or missing opacities) and/or a miscalibration of the SpT-Teff scale. The simplest modification 
to the models (changing to match observations) would yield an older age for this system, in line 
with the recently proposed older age of Upper Scorpius (r ~ 11 Myr). 

Subject headings: 


1. INTRODUCTION 

The fundamental properties of stars constitute a 
bedrock upon which much of astronomy is built, but 
there remains a paucity of well-characterized stars for 
calibrating stellar evolutionary models at young ages (1- 
10 Myr), particularly for low masses (<1 Mq). It is 
a maxim of stellar astrophysics that the properties and 
lifecycle of a star are largely set by its mass, and hence 
it is crucial to calibrate the mass predictions of mod¬ 
els. Uncertainties in model-derived properties constitute 
the dominant systematic error for in-situ measurements 
of the IMF for young populations (Bastian et al. 2010), 
determinations of star cluster ages (e.g., Preibisch et al. 
2001 versus Pecaut et al. 2012), comparison of proto¬ 
planetary disk properties to those of mature planetary 
systems (Andrews et al. 2013), and binary formation 
studies (e.g., Kraus et al. 2011a; Kraus & Hillenbrand 
2012; Duchene & Kraus 2013. Mass/age ambiguities 
strongly limit interpretations for directly-imaged gas gi¬ 
ant planets in young populations, where the host-star 
age and companion luminosity establish the planet mass 
(e.g., Lafreniere et al. 2008 and Ireland et al. 2011 versus 
Pecaut et al. 2012, and also Carson et al. 2013 versus 
Hinkley et al. 2013). 

At present, different model sets predict masses that dif¬ 
fer by 50% for nominally solar-mass pre-main-sequence 
stars at the same point in the L-Tgff HR diagram (e.g., 
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Hillenbrand & White 2004; Torres et al. 2010; Gennaro 
et al. 2012; Stassun et al. 2014a). Furthermore, even the 
temperature scale for converting spectral type to Teff is 
uncertain for young stars at the level of 100-200 K (Luh- 
man et al. 2003; hereafter L03), constituting a systematic 
uncertainty on HR diagram positions. Stellar evolution¬ 
ary models can be tested indirectly, such as by placing 
binary systems (White et al. 1999; Kraus & Hillenbrand 
2009; Torres et al. 2013) or even entire stellar populations 
(Preibisch et al. 2002; Naylor 2009; Pecaut & Mamajek 
2013; Chen et al. 2014) on the HR diagram. However, 
these tests are agnostic to stellar mass and age; they 
demonstrate whether there are any common ages and 
masses that reproduce the HR diagram positions, but 
not whether those values are actually correct. Dynami¬ 
cal masses for young visual binaries have cast much light 
on this issue (Boden et al. 2005; Schaefer et al. 2012; Si¬ 
mon et al. 2013), and are starting to become available in 
statistically robust samples Schaefer et al. 2014; Rizzuto 
et ah, submitted). However, radius measurements are 
more time-dependent (particularly among M dwarfs that 
largely evolve down the Hayashi track to the ZAMS) and 
have remained elusive. Only a small number of young, 
low-mass eclipsing binaries are known (van Eyken et al. 
2011; Morales-Calderon et al. 2012), and with a few well- 
characterized exceptions (Stassun et al. 2006; Irwin et al. 
2007; Cargile et al. 2008; Gomez Maqueo Chew et al. 
2012; Gillen et al. 2014), those systems are faint and re¬ 
main only coarsely characterized. 

The ages of stars remain similarly uncertain, and are 
even more difficult to measure without appeals to evolu¬ 
tionary models. Traceback simulations have been noto¬ 
riously contentious (Ortega et al. 2002; Makarov et al. 
2005; Mamajek & Bell 2014), and HR diagram ages 
can differ significantly depending on the stellar mass 
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Fig. 1. — Left: Postage stamp of the young star UScoCTIO 5 that was downloaded as part of K2’s Campaign 2, rotated so that 
north is up. The red circle shows the 1.5 pixel photometric aperture used in our analysis, after centroiding. The green dot marks the 
nominal location of the star based on target pixel file header position information. Right: A DSS R-band postage stamp of UScoCTIO 5 
(FOV=60'^ North=Up) showing the K2 postage stamp (red dotted box) and our adopted photometric aperture (red circle). The image is 
shown in a square-root stretch using the CubeHelix color palette (Green 2011). All nearby sources fall outside our photometric aperture. 
We have previously shown from analysis of 2MASS images that there are no detected sources with a contrast of AKs < 3 at p > 1.5" and 
AKs < 4 at p > 2.0"(Kraus &: Hillenbrand 2007). UScoCTIO 5 has not been observed with adaptive optics, so there are no limits on 
closer companions, though the absence of a third set of spectral lines in our spectra suggests that there are no additional objects within 
AR <1-2. 


range considered and the assumed star formation history 
(Preibisch et al. 2002; Pecaut et al. 2012; Rizzuto et al. 
2015). The ongoing debate over the age of the Upper 
Scorpius OB association and other young populations 
(Naylor 2009; Bell et al. 2013; Mamajek & Bell 2014) 
demonstrates that systematic errors remain at the level 
of factors of 2. However, if an evolutionary model re¬ 
produces all of the other observable parameters of a star 
(mass, radius, Teff, and Lhoz), then the model-derived 
age could represent the most robust possible theoretical 
estimate of the age for that stellar population. 

To address these fundamental issues of stellar astro¬ 
physics, and to look for planets, we have initiated a 
search for eclipsing/transiting systems among all of the 
known and suspected members of Upper Scorpius and 
Ophiuchus that fell in the footprint observed by K2, the 
extended Kepler mission (Howell et al. 2044) during its 
Campaign 2. These observations comprised 79 contin¬ 
uous days of observing with a 30 minute cadence, with 
no gaps, and hence are ideal for identifying all eclipsing 
binaries with periods on the order of this duration or 
shorter. In total, we proposed observations of 657 con¬ 
firmed members (GO2052, PI Covey) and 759 likely can¬ 
didate members (GO2063, PI Kraus) that were optically 
visible {Kp < 16 mag). Given the size of our sample and 
the known frequency of short-period binary systems in 
the field (Raghavan et al. 2010), we expect that ^10-15 
new eclipsing binaries should be discovered. 


The young Upper Scorpius member UScoGTIO 5 (Fig¬ 
ure 1)® was a high-priority target in our program since it 
is a known spectroscopic binary, first reported and char¬ 
acterized as such by Reiners et al. (2005). The member¬ 
ship of UScoGTIO 5 in Upper Scorpius was first proposed 
by Ardila et al. (2000), based on its HR diagram position, 
line of sight extinction, and the presence of Ha emis¬ 
sion; this membership was confirmed by Reiners et al. 
(2005) based on the presence of lithium absorption that 
constrains the age to be r < 30 Myr. Reiners et al. 
(2005) also detected the presence of two sets of spectral 
lines in their first high-resolution spectrum, and subse¬ 
quently obtained enough additional spectroscopic obser¬ 
vations to determine the period {P = 34 days), the or¬ 
bital elements, and a minimum mass (modulo the sine 
of the orbital inclination). Their inferred minimum mass 
{Mtot sin(*) = 0.64 ± 0.02 Mq) was larger than the mass 
predicted by evolutionary models, which suggested both 
that the models could require modifications or additional 
physics, and that the system could be close to edge-on 
(and hence could show eclipses). 

In this paper, we report that UScoCTIO 5 is indeed 
oriented to show both primary and secondary eclipses in 

® The USco candidates identified by Ardila et al. (2000) are 
sometimes abbreviated as “USco NN”. However, SIMBAD has 
assigned the name “USco 5” to different candidate member by 
Sciortino et al. (1998). We caution the reader to not confuse the 
two sources if pursuing additional followup observations. 
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its K2 light curve. We re-interpret the sum of data avail¬ 
able for the system to determine masses and radii for each 
star, to test different evolutionary model tracks against 
the empirical constraints of this system, to calculate an 
empirical constraint on the SpT-Teff conversion for young 
mid-M stars, and to determine a new semi-empirical age 
estimate for the Upper Scorpius OB association. 

2. K2 LIGHT CURVES 
2 .1. K2 Photometry 

We downloaded the target pixel hies (TPFs) for US- 
coCTIO 5 from the Mikulski Archive for Space Tele¬ 
scopes (MAST), where it is stored under its Eclip¬ 
tic Plane Input Catalog (EPIC) identihcation number 
205030103. This data consists of 3811 ~10xl2 pixel 
stamp images centered on UScoCTIO 5, acquired be¬ 
tween 2014 August 23 and 10 November. 

Telescope pointing for K2 is unstable, due to the fail¬ 
ure of a second reaction wheel during the Kepler prime 
mission in May of 2013. As a result, the centroid of 
UScoCTIO 5 drifts at a rate of ^0.1", or ~0.02 pixels, 
per hour. While this movement is relatively small, asso¬ 
ciated detector sensitivity variations at the few percent 
per pixel level compromise the otherwise exquisite pho¬ 
tometry. To compensate for these effects, onboard soft¬ 
ware checks the pointing every six hours and initiates a 
thruster bring if the roll angle has exceeded a threshold. 
The result is jumps in position (and hence measured hux) 
on this timescale. 

To minimize the effect of jumps on the photometry, we 
determined the stellar centroid position independently 
for each pixel stamp. TPF headers provide a rough world 
coordinate system solution which is the same for all im¬ 
ages, but these are not precise enough to center the tar¬ 
get. We therefore cut out a 7x7 pixel substamp around 
the nominal target position (i.e., the green dot in the left 
side of Figure 1), and determined a flux-weighted cen¬ 
troid from these pixels. The drift in these centroids over 
time can be seen in the bottom panels of Figure 2. 

We then placed photometric apertures at each centroid 
location, extracting photometry with a set of circular 
apertures from 1.0 to 4.0 pixels in radius, at intervals of 
0.5 pixels, and subtracting the background as determined 
from a wider annulus. The circular aperture is assumed 
to have a boundary that intersects the square pixels along 
an arc, with fractional hux per pixel integrated geometri¬ 
cally for pixels neither fully inside nor outside the circu¬ 
lar aperture. We found that photometric noise levels af¬ 
ter detrending for position jump effects (see Section 2.2) 
were minimized with the 1.5-pixel aperture. This size 
has the additional advantage of being small enough so as 
to avoid hux contamination from other stars lying ^10" 
away. The chosen aperture is shown as a red circle in 
Figure 1. 

To remove errant data, we discarded the hrst 93 light 
curve points, for which the pixel positions were particu¬ 
larly different compared to the rest of the time series. We 
also removed points with detector anomaly hags. Finally, 
we pruned points lying more than hve standard devia¬ 
tions off the median light curve trend (excluding points 
within or around eclipses). The resulting raw light curve 
is displayed in Figure 2. 


2 .2. Detrending of Instrumental Effects and Stellar 
Variability 

Before eclipse htting was performed, the K2 light curve 
of UScoCTIO 5 was corrected for instrumental artifacts 
and detrended to remove starspot signatures that domi¬ 
nate UScoCTIO 5’s out-of-eclipse light curve. The degra¬ 
dation of Kepler’s pointing control in the K2 mission re¬ 
sults in light curves which, in many cases, are dominated 
by the target star’s path across small (1-2%) sensitiv¬ 
ity variations in Kepler’s detector. However, as Van- 
derburg & Johnson (2014) and Vanderburg et al. (2015) 
have shown, the well-behaved nature of this correlation 
allows the removal of much of this signature by applying 
a position-dependent correction to each K2 light curve. 

We applied such a correction to UScoCTIO 5’s light 
curve, removing the systematic structure associated with 
the drift of the K2 focal plane. To infer the correlation 
between UScoCTIO 5’s xy position on the detector and 
the sensitivity of the K2 light curve, we hrst performed 
an initial detrending of UScoCTIO 5’s light curve to re¬ 
move the intrinsic stellar variations that is ubiquitous 
to young stars (e.g., Cody & Hillenbrand 2010; Cody 
et al. 2013, 2014) and whose amplitude is lOx larger 
than the systematic pointing artifacts. We removed 
these stellar variations using the SuperSmoother algo¬ 
rithm, originally developed by Friedman (1984) and sub¬ 
sequently implemented in python by Vanderplas (2015; 
10.5281/zenodo. 14475), by dividing the raw light curve 
by an alpha=0 supersmoothed ht. The alpha parameter 
provides a mechanism for biasing the SuperSmoother ht 
against high frequency components (i.e., providing ’bass 
enhancement’ in Friedman’s original treatment); an al- 
pha=0 ht corresponds to an unbiased ht, which preserves 
the ability to ht and correct for high frequency stellar 
variations in the target light curve. 

With the large-scale stellar variability signatures re¬ 
moved, the systematic structure due to K2’s pointing 
drift dominated the normalized light curve; we removed 
this structure by dividing the hux in each epoch of US¬ 
coCTIO 5’s raw light curve by the median normalized 
hux of the 10 data points in the normalized light curve 
whose xy positions are closest to that of the epoch in 
question. The correlation between K2’s sensitivity and 
UScoCTIO 5’s xy position underwent a clear change 
halfway through campaign 2, when the direction of the 
torque resulting from the incident solar hux changed to 
the opposite direction. As a result, we corrected the por¬ 
tions of the light curve taken before and after epoch 2102 
independently.^ 

Finally, having applied this pointing correction to USco 
CTIO 5’s raw light curve, we then again used an a = 0 
SuperSmoother ht to detrend the light curve and remove 
the large-scale starspot signatures prior to htting the 
eclipse prohles. The remaining light curve still showed 
long-term variations at the level of ±3 millimagnitudes, 
so we removed this remaining small signal by htting (us¬ 
ing datapoints outside eclipse) and subtracting a DC off¬ 
set for ±12 hours around each eclipse. We show the 
raw, instrumental-detrended, and instrumental/stellar- 
detrended light curves in Figure 3. 

^ In the text and in Figures 2 and 3, we use the Kepler time 
coordinate system for the K2 light curve where epoch = BJD - 
2454833. 
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Fig. 2.— Aperture photometry results for the young star UScoCTIO 5. Time is specified in units of BJD-2454833. Top: Normalized 
light curve extracted from aperture photometry, without any subsequent detrending. The eclipse events can be seen at days 2077, 2088, 
2110, and 2121. Middle and Bottom: X and Y centroid positions, in pixels, as a function of time. The six-hour interval between thruster 
firings (which reset the telescope position) is evident in the positions, and the position information can be used to detrend flux variations 
as the target moves across the detector. 


3. HIGH-RESOLUTION SPECTROSCOPY AND RADIAL 
VELOCITIES 

High-resolution spectroscopic measurements for US- 
coCTIO 5 were already reported by Reiners et al. (2005), 
who obtained 22 epochs in 2003-2005 with the optical 
echelle spectrographs Keck-I/HIRES, Magellan/MIKE, 
and Blanco/Echelle. They extracted the difference be¬ 
tween the primary and secondary star RVs (Av = Vp—Vg) 
for the system in each epoch, determined from a cross¬ 
correlation with the dwarf standard star G1 406. They 
reported a best-fit orbit with P = 33.992 ± 0.006 days, 
(Ml -P M 2 )sin(t) = 0.64 ± 0.02Mq, e = 0.276 ± 0.008, 
a = 0.177 ± 0.002 AU, w = 355.5 ± 0.8 deg, and 
To = 52799.974 ± 0.002 days (MJD). However, since the 
K2 light curve offers precise constraints on the period 
and a combination of the eccentricity and longitude of 
periastron, we have reht the orbit as part of our analy¬ 
sis. 

We have adopted the Magellan and Blanco Av mea¬ 
surements directly from the analysis by Reiners et al. 
(2005). Based on the RMS scatter that they reported for 
their measurements around the best-fit orbit, the typical 
uncertainty in each measurement is (tau = 500 m/s. Eor 
the Keck-I/HIRES data, we have downloaded the obser¬ 
vations from the Keck archive and reanalyzed the spectra 
to measure absolute RVs for each component. These ab¬ 
solute velocities are needed in order to measure the mass 
ratio of the system, and hence the masses of the individ¬ 
ual component stars. 


Our analysis of the HIRES data is very similar to the 
methods described in Kraus et al. (2011b) and Kraus 
et al. (2014). We extracted and wavelength-calibrated 
each spectrum using the MAKEE pipeline®, refining 
the wavelength solution by cross-correlating the 6800A 
or (preferably) 7600A telluric absorption bands against 
those of the spectrophotometric standard star HZ 44 
(Massey et al. 1988). Eor each spectrum of UScoCTIO 
5, we then measured the broadening function (Rucinski 
1999)® with respect to our own Keck/HIRES observa¬ 
tions of two standard-stars over a total of three epochs: 
one observation of G1 447, and two separate observa¬ 
tions of G1 83.1. Both standard stars have temperature 
and metallicity similar to UScoCTIO 5 (SpT=M4-M5 
and [M/H] ^ 0.0; Rojas-Ayala et al. 2012, Mann et al. 
2015b). The broadening function is a better represen¬ 
tation of the rotational broadening convolution than a 
cross-correlation, since it is less subject to “peak pulling” 
and produces a flatter continuum. We ht each broaden¬ 
ing function with two Gaussian functions to determine 
the absolute primary and secondary star RVs (vp and 
Vg), the standard deviations of the lines due to rotation 
and instrumental resolution (ct^^ and a-y ^), and the aver¬ 
age flux ratio across all well-fit orders (which is estimated 
from the ratio of areas for the two peaks of the broaden¬ 
ing function). 

We list these measurements in Table 1. Eor those spec- 

® http://spider.ipac.caltech.edu/stafF/tab/makee 

® http://www.astro.utoronto.ca/'rucinski/SVDcookbook.html 
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Fig. 3.— Light curve for UScoCTIO 5 against time in units of BJD-2454833. Top: The raw fluxes measured via aperture photometry. 
Middle: Corrected fluxes where the detector-position-dependent brightness changes have been detrended out. Bottom: Corrected fluxes 
where the out-of-eclipse variations due to spots have been detrended out. The variability is clearly dominated by spot-driven variations on 
a characteristic timescale of ~34 days and with amplitude of ±5%, albeit with clear changes in shape between subsequent cycles. The two 
primary and two secondary eclipse events, which have depths of ~5% and 7% respectively, are clearly seen. Based on detailed inspection 
of the light curve, the numerous brightening events appear to be astrophysical (stellar flares) and not systematic errors. 


TABLE 1 

Keck-I/HIRES RVs 


Epoch 

Epoch 

Wavelength 

^int ) 

SNR^ 

Vp 


Vs 


Fs/Fp 

(UT Date) 

(HJD-2450000) 

Range (A) 

(sec) 


(km/s) 

(km/s) 

(km/s) 

(km/s) 


20030611 

2801.94316 

6400-8700 

900 

49 

-36.78± 0.18 

8.05± 0.10 

33.22± 0.17 

8.09± 0.27 

0.915± 0.034 

20040509 

3134.96680 

7000-9250 

900 

50 

-9.18± 0.12 

6.65± 0.32 

4.29± 0.13 

7.18± 0.24 

1.045± 0.052 

20040510 

3136.00935 

7000-9250 

1200 

52 

-15.49± 0.09 

6.92± 0.22 

11.47± 0.21 

7.01± 0.27 

0.941± 0.011 

20040511 

3136.90724 

6400-8700 

1200 

43 

-22.05± 0.16 

7.71± 0.23 

17.25± 0.19 

7.27± 0.28 

0.908± 0.026 

20040511 

3137.00804 

6440-8750 

900 

39 

-22.53± 0.10 

7.80± 0.14 

18.27± 0.11 

7.66± 0.14 

0.932± 0.020 

20040511 

3137.13712 

6440-8750 

400 

5 

-24.44± 0.17 

7.74± 0.21 

18.57± 0.16 

7.52± 0.31 

0.951± 0.044 

20040513 

3139.13254 

4480-6890 

400 

17 

-35.89± 0.20 

7.59± 0.21 

31.65± 0.20 

8.04± 0.19 

0.932± 0.034 

20050302 

3432.06727 

5720-8570 

900 

33 

17.79± 0.13 

7.85± 0.13 

-24.12± 0.11 

7.87± 0.12 

0.922± 0.017 

G1 83.1 

5741.13595 

4320-8750 

120 

46 






G1 83.1 

5930.69346 

4320-8750 

120 

71 






G1 447 

5933.16968 

4320-8750 

120 

130 






HZ 44 

5931.18067 

4320-8750 

120 

43 







^Measurements of the spectrum’s SNR and the components’ flux ratio are made at A = 6600A if that wavelength is included in the 
observation, and at A = 7000A otherwise. 

^We report (Jvp and cry^ as the standard deviation of the Gaussian fits to the two stars’ broadening functions, which is a measure of 
both the instrumental broadening and the rotational broadening. 


tra that include the 6500-6700A region, we also report 
equivalent widths of Ha emission (which was never fully 
resolved due to the intrinsic thermal broadening of the 
chromosphere, and hence we report as a single value) and 
Ligyos absorption (which was resolved in all epochs) in 
Table 2. The equivalent widths are measured with re¬ 


spect to the continuum of the full composite spectrum, 
but individual stellar values can be determined from the 
flux ratio of the spectra (which is nearly constant across 
the entire wavelength range of these spectra). We also 
show narrow wavelength ranges around Ha and Ligyog in 
each spectrum in Figure 4. 
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Fig. 4.— Hq emission and lithium absorption lines from the six Keck/HIRES spectra that included the appropriate wavelength ranges. 
All spectra are normalized to the (pseudo-)continuum surrounding that line. All were observed near quadrature, resulting in clearly resolved 
lithium lines. However, the Ha emission lines are never more than moderately resolved; they are broader because they are formed in the 
hot chromosphere, unlike the photospheric absorption lines that are only rotationally broadened. 


TABLE 2 

Keck-I/HIRES Spectral Features 


Epoch 
(UT Date) 

HJD 

(days) 

EW[Ha] 

(A) 

E W [ _L 2 ] p 7. j 772 

(A) 

EW[Li]s<,c 

(A) 

20030611 

2801.94316 

-6.1 

0.32 

0.23 

20040511 

3136.90724 

-4.8 

0.29 

0.27 

20040511 

3137.00804 

-3.8 

0.31 

0.29 

20040511 

3137.13712 

-4.0 

0.25 

0.20 

20040513 

3139.13254 

-4.7 

0.32 

0.27 

20050302 

3432.06727 

-4.6 

0.34 

0.30 


4. SPECTRAL CLASSIFICATION, EXTINCTION, 
TEMPERATURE, AND LUMINOSITY 

4.1. Intermediate-Resolution Speetroscopy 

To better determine its spectral type and extinction, 
we observed UScoCTIO 5 on 2015 April 1 using the 
Wide-Field Spectrograph (WIPES) on the ANU 2.3ni 
telescope. These observations are identical to those de¬ 
scribed by Rizzuto et al. (2015), which describes the data 
and procedures in more detail. WIPES is an integral field 
spectrograph with an FOV of 38x25", feeding red and 
blue arms. We configured the red arm to deliver a spec- 
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Fig. 5.— Left: The WIPES optical spectrum of UScoCTIO 5 as compared to a range of dwarf spectral standards. In each case, we fit 
for the extinction Ay that minimizes the residuals; the best-fit spectrum is an M4.5V dwarf with Ay = 1.1 added. Right: The x? surface 
describing an expanded range of fits like those shown. The blue dots show the best-fit value of Ay for each spectral type. The contours 
are drawn at at levels of Axu =1, 2, 3, 4, and 5; it is difficult to infer accurate uncertainties from fits to spectra due to the strong 
covariances in both astrophysical and instrumental errors, but the Ax? = 1 contour demonstrates the range over which the models visually 
give a worse fit. From this and from visual inspection of the left panel, the uncertainties are ±0.5 subclass in SpT and ±0.3 mag in Ay. 


TABLE 3 

System Photometry 


Filter 

m (mag) 

Reference 

B 

17.806 ± 0.203 

APASS (Henden et al. 2012) 

V 

16.192 ± 0.010 

APASS (Henden et al. 2012) 

g' 

16.975 ± 0.074 

APASS (Henden et al. 2012) 

r' 

15.482 ± 0.042 

APASS (Henden et al. 2012) 

r' 

15.385 ± 0.050 

CMC15 (Evans et al. 2002) 

i' 

13.708 ± 0.011 

APASS (Henden et al. 2012) 

J 

11.172 ± 0.023 

2MASS (Cutri et al. 2003) 

H 

10.445 ± 0.026 

2MASS (Cutri et al. 2003) 

Ks 

10.170 ± 0.021 

2MASS (Cutri et al. 2003) 

W1 

10.036 ± 0.023 

ALLWISE (Cutri & et al. 2013) 

W2 

9.838 ± 0.020 

ALLWISE (Cutri & et al. 2013) 

W3 

9.648 ± 0.047 

ALLWISE (Cutri & et al. 2013) 


tral resolution of R = 7000 across a wavelength range of 
A =5300-7000A; the flux to the blue arm was not sufh- 
cient to deliver useful data. We observed UScoCTIO 5 
with an integration time of t = 500 s, yielding S/N = 50 
at A =6600A. We processed the raw WIPES data with 
the “WiFeS PyPeline” (Childress et al. 2014) software 
in order to extract a spectral data cube. We then mea¬ 
sured the flux in each spectral channel using PSF-fitting 
photometry with a Moffat profile, integrating the pro¬ 
file to measure the total flux from the science target in 
each channel and subtracting the sky background implied 
from the DC offset of the PSF fit. We show the extracted 
spectrum in Figure 5 (left); for further detail regarding 
the extraction process see Rizzuto et al. (2015). 
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Fig. 6. — Photometric SED of the (unresolved) UScoCTIO 5 
system (red points), as compared to the reddened best-fitting BT- 
SETTL model (T = 3250 K, Ay =0.9) and the corresponding syn¬ 
thetic photometry in those bands. The bottom panel shows the 
difference between the synthetic and observed photometry in stan¬ 
dard deviations. 


4.2. Archival Photometry 

The geometric determination of radii offers the intrigu¬ 
ing possibility of directly measuring empirical tempera¬ 
tures, as long as an accurate luminosity can be deter¬ 
mined. While it is possible to estimate the stellar lumi¬ 
nosity from a single flux and a corresponding bolomet- 
ric correction, it should be more accurate to compile a 
broadband SED and add up the flux. This strategy also 
minimizes model-dependent uncertainties in the assumed 
























bolometric correction. To that end, we also have com¬ 
piled all of the available (component unresolved) pho¬ 
tometry in all-sky surveys. As we summarize in Table 3, 
we have used photometry from 2MASS {J,H,Ks; Cutri 
et al. 2003), AllWISE {Wl, W2, W3; Cutri & et al. 
2013), CMC15 (r'; Evans et al. 2002), and APASS {B, 
V, g', r', i'; Henden et al. 2012). 

4.3. Atmospheric Properties 

We have used the low-resolution optical spectrum from 
WIPES (Section 4.1; Eigure 5, left) to calculate a joint 
constraint on the spectral type and extinction of the US- 
coCTIO 5 system. We compared the observed spectrum 
to a sequence of field M dwarfs compiled from SDSS 
(Bochanski et al. 2007), that we artificially reddened us¬ 
ing an Ry = 3.1 reddening law (Savage & Mathis 1979); 
the best spectrum is an M4.5 dwarf that has been red¬ 
dened by Ay = 1.06 mag. In Eigure 5 (right), we show 
the best-fitting Ay value as a function of spectral type 
for the full grid of standards, as well as contours in the 
xt surface. It is not straightforward to convert these con¬ 
tours to a statistically robust confidence interval because 
spectra show strong covariances across wide wavelength 
ranges due to both astrophysical and instrumental er¬ 
rors. However, the interval of = 1 denotes the range 
across which the residuals become visually apparent, and 
hence represent an upper limit on the uncertainty in the 
spectral classification. Since spectral types themselves 
are only defined by half-subclasses, then we assess a final 
uncertainty of ±0.5 subclass on the spectral type and a 
corresponding uncertainty of lto’27 mag on the extinction. 

The conversion from spectral type to temperature has 
been a matter of longstanding debate in the star forma¬ 
tion community. It is well known that this conversion is 
gravity-sensitive for M dwarfs (e.g., Luhman 1999) and 
perhaps even for earlier-type stars (Pecaut & Mamajek 
2013), such that giants of equivalent spectral type are 
several hundred degrees hotter than dwarfs. L03 pro¬ 
posed an “intermediate gravity” temperature scale that 
might be appropriate for young stars, constructed from 
an interpolation of the dwarf and giant scales that makes 
the GG Tau system look coeval (as per White et al. 
1999). Given the spectral type of M4.5±0.5, the in¬ 
ferred temperature on this intermediate-gravity scale is 
Teff = 3200 ±75 K. The dwarf scales would imply a lower 
temperature of Teff = 3100 ± 90 K according to Leggett 
et al. (1996), or warmer temperatures of Tes = 3190±75 
K from more recent scales (Rajpurohit et al. 2013; Mann 
et al. 2015b) while the giant scale (van Belle et al. 1999) 
would imply a higher temperature of Teff = 3435 ± 55 
K. Herczeg & Hillenbrand (2014) more fully review the 
wide range of other proposed temperature scales, which 
yield temperature estimates spanning 3080-3305 K for 
an M4.5 dwarf star; they develop an independent young 
star temperature scale that would predict Teff = 3085 K 
for UScoCTIO 5. Eor this work, we adopt the temper¬ 
ature scale of L03 and the corresponding spectroscopic 
temperature of Teff = 3200 ±75 K in order to remain con¬ 
sistent with the previous literature, but we also remind 
the reader that substantial systematic uncertainties re¬ 
main. 

As we discuss further in Section 6.2, it is possible to em¬ 
pirically anchor the SpT-Teff relation at M4.5 using our 
geometric radius measurements if we can also measure 


the luminosity. To that end, we computed Ftoi follow¬ 
ing the procedure from Mann et al. (2015b). We used a 
set of solar-metallicity, log((7) = 4.25 BT-SETTL models 
(Allard et al. 2011) built with the Caffau et al. (2011) 
solar abundances, which better reproduce the observed 
spectra of Held M dwarfs than BT-SETTL models utiliz¬ 
ing the Asplund et al. (2009) abundances (Mann et al. 
2013). We selected grid points spanning the range of 
temperatures found above, each of which we reddened 
by a grid of extinction values consistent with the value 
and error derived using the empirical spectrum. We then 
scaled each model spectrum to give the best agreement 
(lowest x^) between synthetic and observed photometry 
and integrate over the (un-reddened but scaled) model 
spectra to determine Tho/- The range of temperatures 
and extinction values produced a range of Fhoi values, 
which we use as an estimate of the uncertainty. The 
best fit model {Tmodei = 3250 K, Ay = 0.9 mag) is 
shown in Figure 6. To convert from Fboi to Lboi, we 
assume a distance of d = 145 ± 15 pc, consistent with 
the mean and scatter seen for BAF stars in Upper Scor- 
pius (e.g., de Zeeuw et al. 1999; Rizzuto et al. 2011). 
The ~20% uncertainty in extinction dominates the un¬ 
certainty in Fboi {fJFboi ~5%). The 10% uncertainty in 
the distance, which propagates to a 20% uncertainty in 
the luminosity, represents the dominant source of un¬ 
certainty in Lboi- Using the model to fill the gaps be¬ 
tween the photometry, we find a total bolometric flux of 
Fboi = (2.02tQ;|!)g) x 10“^° erg/s/cm^. For an assumed 
distance of d = 145 ± 15 pc, the corresponding lumi¬ 
nosity is Ltot = O-132j'jQ g3gT0, or individual component 
luminosities of T ~ O.O65T0. 

5. RESULTS 

5.1. Fitting for Orbital and Stellar Parameters 

In the simple case of a tidally circularized orbit, fitting 
for the parameters of an eclipsing binary is a straight¬ 
forward and separable problem: the component masses 
(moduli sin(i)) can be derived from the radial velocity 
curve, while the inclination, period, radii, and stellar flux 
ratio can be derived from the light curve. However, ec¬ 
centric systems are more complicated, with many covari¬ 
ances between parameters. For example, while the RV 
curve provides a direct measurement on the eccentricity 
e and the argument of periapse w, the orbital phases of 
primary and secondary eclipse can provide a very tight 
joint constraint on both parameters. The durations of 
the eclipses provide a further constraint on e and w, as 
well as posing a joint constraint on the inclination and 
surface brightness ratio (since the surface areas occulted 
can differ between the primary and secondary eclipses). 
Finally, since the RVs and photometry were measured 
nearly a decade apart, then timekeeping errors could lead 
to covariances in the parameters. For a system this wide, 
apsidal motion should be negligible on decade timescales 
(e.g., Feiden & Dotter 2013), so we do not expect the 
orbital parameters to have changed between the epochs 
of the RV data and K2 data. 

To properly account for these covariances, we have con¬ 
structed a Markov-Chain Monte Garlo (MGMC) proce¬ 
dure to simultaneously fit both the RV curve and the 
light curve using a full model of all parameters. We 
specifically consider a model where the explicitly fit pa¬ 
rameters are 6 orbital elements (Tq, P, a, e, w, and i, 
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Fig. 7. — Left: Absolute radial velocities Vp and Vs for the primary and secondary stars of UScoCTIO 5, as measured from the 
Keck/HIRES epochs listed in Table 1. We also show the best-fit model as determined from our fitting procedure (Section 5.1). Underneath, 
we show the (0-E) residuals with respect to the best-fit model. Right: Relative velocity differences Av = Vs — Vp for all epochs, including 
the measurements reported by (Reiners et al. 2005), also with the best-fit model curve and (O-E) residuals shown. 
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Fig. 8 .— Primary eclipse (left) and secondary eclipse (right) for UScoCTIO 5, along with the best-fitting models (dashed lines) and 
the (O-E) residuals (bottom panels). Data points that were rejected as outliers are shown with open circles; a flare occurred during one 
secondary eclipse, so we have rejected all affected points. The primary eclipse is deeper primary due to the orbital geometry, such that a 
larger total surface area is occulted when the primary star is behind the secondary star; the surface brightnesses of the two stars are nearly 
identical. Despite the grazing nature of the eclipses, the full radii can be determined with small uncertainties and minimal covariance due 
to constraints from the timing of primary versus secondary eclipse, the relative eclipse durations, the relative eclipse depths, and the precise 
RV curve. 


neglecting fl in the absence of spatially resolved infor¬ 
mation), the mass ratio of the system q = Mg/Mp, the 
systemic radial velocity 7 , the sum of the stellar radii 
Rtot = Rp + Rg, the ratio of the stellar radii r = Rg/Rp, 
and the ratio of stellar fluxes through the Kepler Kp 


bandpass f = Fg/Fp. 

We note that the convention for eclipsing binaries is 
not to fit for / = Fg/Fp, but rather for a temperature 
ratio or temperature difference (indeed, as we have done 
in the past; Kraus et al. 2011b). In principle, the ratio 
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TABLE 4 


System Parameters for UScoCTIO 5 

Orbital Parameters 

To (HJD) 

P (days) 
a (AU) 
e 

i (deg) 
uj (deg) 

7 (km/s) 

56914.490 ± 0.026 
34.00073 ± 0.00007 
0.17749 ± 0.00031 
0.26741 ± 0.00011 
87.912 ± 0.010 
355.13 ± 0.30 
-2.64 ± 0.07 

Stellar Bulk Parameters 

Mp + Ms (Mo) 
q = Ms /Mp 

Mp (Mo) 

Ms (Mo) 

Pp “t Ps {Pq) 

Ps/Pp 

Pp (Pq) 

Ps {Pq) 
log(ffp) (cm/s2) 
log(ffs) (cm/s^) 

0.6452 ± 0.0034 
0.963 ± 0.007 
0.3287 ± 0.0024 
0.3165 ± 0.0016 
1.644 ± 0.008 
0.972 ± 0.011 
0.834 ± 0.006 
0.810 ± 0.006 

4.11 ± 0.01 

4.12 ± 0.01 

Stellar Atmospheric 

Parameters 

Fs/Fp 

Ts/Tp 

Av (mag) 

SpT 

Piff.LOS (K) 
Piff.HHlT (K) 

Fbol (erg/s/cm^) 

0.999 ± 0.017 

1.000 ± 0.004 

1 06"*"®'^^ 

27 

M4.5 ± 0.5 
3200 ± 75 
3085 ± 105 
(2.02l®;J|) X 10-1® 

Parameters Using Distance 

Lboi (-^o) 

Lboi (Lq) ' 

^effjgeom (K) 
^eff,geom (K) 

n 1 oo”f0.028 

U.iUZ_Q Q30 

0.132l®;°®%D/145pc)2 

O 9 OK + I 6 O 

ozoo_ 2 oo 

323511° (D/145pc)2 

Infrared f'lux Method Distance 

dL03 (pc) 
dHH14 (pc) 

144.4 ± 6.6 
135.1 ± 8.8 


of surface brightnesses can be described by the ratio of 
temperatures; in the limiting case of a pure blackbody, 
the relation is exact and analytic. However, cool stellar 
atmospheres are notably non-blackbody, making this re¬ 
lation more complex. Traditional fitting codes like the 
Wilson-Devinney algorithm (Wilson & Devinney 1971, 
and extensive updates thereof) parametrize this relation 
using Kurucz model atmospheres (Kurucz 1979). How¬ 
ever, those codes typically only extend to T ~ 3500 K 
(hotter than our observed targets), and the conversion 
from temperature ratio to flux ratio is not easily quan¬ 
tified or changed. We find it more straightforward to fit 
directly for the flux ratio in the Kepler bandpass, and 
then deal with the conversion to a temperature ratio ex¬ 
plicitly. This choice can also be found in other fitting 
codes such as JKTEBOP (Southworth et al. 2004), which 
shares many design choices with our procedures. 

In executing our MCMC procedure, we use analytic 
equations to construct a predicted RV curve against 
which we can compare the observations and measure the 
residuals. In order to isolate potentially correlated un¬ 
certainties between the measurements of the primary and 
secondary star at each epoch (such as from uncertainties 
in the wavelength scale), we chose to fit the primary star 
RV Vp (for the 8 epochs of HIRES data) and the differ¬ 
ence in RV Av = Vs—Vp (for all 22 epochs). Some fit pa¬ 
rameters (such as the semimajor axis, and hence the total 
system mass) only depend on Av, while others (such as 


the system velocity and the mass ratio) necessarily de¬ 
pend on the individual component RVs, so this choice 
ensures that the correlated RV errors between primary 
and secondary star will not unnecessarily inflate the un¬ 
certainty in parameters that do not require component- 
resolved measurements. 

The analysis of the light curve is less straightforward. 
As we mention above, the Wilson-Devinney code is the 
gold standard of the field due to the wide range of phys¬ 
ical effects (such as reflected light and tidal distortion) 
that it can reproduce. However, most of these physi¬ 
cal effects are not needed for a binary with a semima¬ 
jor axis of nearly 0.2 AU, and the architecture needed 
to encompass them results in a long runtime (of order 
seconds) to produce a single model. Eurthermore, the 
Wilson-Devinney solution for stars that fall outside their 
modeled temperature range (interpolation between the 
coolest model and a blackbody) is also not appropriate 
for cool stars, meaning that its implicit conversion from 
/ to a temperature ratio would not be correct. 

We therefore have instead constructed an analytic for¬ 
malism that uses the work of Mandel & Agol (2002) to 
calculate the total light removed from the system due to 
occultation of whichever star is more distant. This ana¬ 
lytic model can be calculated ~100 times more quickly, 
allowing for more and longer MCMC chains that better 
explore the complex multi-dimensional parameter space 
of the fit. To account for the long duration of individual 
K2 exposures (which can result in a significant deviation 
between the midpoint flux and the average flux), we cal¬ 
culated the occulted flux at one minute intervals within 
each integration, and then calculated the average value 
for comparison to the observations. We chose our fit pa¬ 
rameters to encompass known covariances (for example, 
fitting Rtot = Rs+Rp and r = Rs/Rp, since the two radii 
are known to be degenerate and anti-correlated in EB 
fits). To allow for limb darkening, we use a quadratic re¬ 
lation with the coefficients prescribed for a star of appro¬ 
priate Teff and log((/) by Claret et al. (2012): 71 = 0.5125 
and 72 = 0.2533. 

Einally, our fit also includes the observed optical flux 
ratio from the Keck/HIRES spectra (0.943 ± 0.015), as 
determined from the ratio of the integrated areas under 
each star’s broadening function. Since the HIRES spec¬ 
tra represent the same wavelength range as the K2 light 
curve, then this comparison directly offers a joint con¬ 
straint on the ratio of surface brightnesses (/ = Fg/Fp) 
and the ratio of radii (r = Rg/Rp). As we mention above, 
the latter parameter in particular can be highly degen¬ 
erate for grazing-eclipse systems without an extra con¬ 
straint. 

We executed the MCMC using a Metropolis-Hastings 
sampler to walk through parameter space, using jump 
sizes drawn from Gaussian distributions with standard 
deviation corresponding to a characteristic jump size. 
We executed several test chains early in this process, 
tweaking the jump sizes to yield acceptance rates of 
20 < P < 50%. We then computed four simultaneous 
chains for a total length of 10® steps per chain. As a re¬ 
sult, our distributions have 4 x 10® distinct samples from 
which the posteriors on each parameter are constructed. 
We also verified that the individual chains yield values 
that agree to within much less than the reported Icr un¬ 
certainties, indicating that they are well-mixed. 
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Finally, we calculated other parameters of interest 
{Mp, Ms, Mtot = Mp + Ms, Rp, Rs, Tes,s/TeB,p) from 
the fit parameters at each step in the chain, yielding 
similar posterior distributions. This method naturally 
propagated the uncertainties and covariances in the fit 
parameters through to the uncertainties in the derived 
parameters. In particular, the ratio of effective tempera¬ 
tures was calculated simply assuming a blackbody, since 
the surface brightness ratio is statistically consistent with 
unity and hence the stellar photospheres (and resulting 
emergent spectra) must be very nearly identical. 

The Mandel & Agol (2002) formalism faces limitations 
for some systems, since it can’t encompass reflected light, 
tidal distortion, or various other effects. However, most 
of these effects should not be relevant for a wide system 
like UScoCTIO 5. The one exception is the influence of 
spots; the 5-10% variations in out-of-eclipse flux suggest 
that both stars are likely to be covered with large and 
complicated spot patterns. If spots are occulted during 
an eclipse, then the different surface brightness of the 
covered area will lead to a different change in the to¬ 
tal system brightness, distorting the eclipse morphology. 
Indeed, we appear to see these effects in our own light 
curve fits at the ±2 mmag level. 

The traditional solution has been to fit with a spot 
model (typically consisting of a few large spots) and op¬ 
timize their size, latitude, longitude, and temperature to 
match the out-of-eclipse brightness. However, these spot 
models are highly degenerate, with many possible con¬ 
figurations replicating the same broad variations. If the 
incorrect spot model is used (as is almost certainly the 
case), then it will degrade the precision of the eclipse fit 
by simultaneously not encompassing the fine details of 
the spot structure (which can not be fit from the varia¬ 
tions in total system flux) and forcing the fit to account 
for a spot model that is not correct. We therefore argue 
that the most conservative solution is to fit with no spots, 
and then forward-model a range of spot models into the 
observational space and determine the resulting scatter 
in best-fit solutions. As we discussed in Kraus et al. 
(20IIb), using this procedure for field M-dwarf eclipsing 
binaries with similar variations resulted in radius uncer¬ 
tainties of ±2%, which we adopt as a systematic uncer¬ 
tainty on our radii in this paper. 

5.2. System Properties 

UScoCTIO 5 is one of the few young low-mass bina- 
ries (r < 10 Myr; O.I < M* < 0.7Mq) for which precise 
masses and radii have been determined, and therefore 
it represents a strong test of pre-main-sequence stellar 
evolutionary models. We summarize our best-fit proper¬ 
ties of UScoCTIO 5 and its component stars in Table 4, 
and in Figures 7 and 8, we show the observed RVs and 
photometry, the best-fit model RV curve and light curve, 
and the residuals between the observations and the data. 
We find that UScoCTIO 5 consists of two nearly iden¬ 
tical components with masses M ~ 0.32Mq and radii 
R ~ O.82i?0; the system mass that we calculate is con¬ 
sistent with the value of M sin(i) measured by Reiners 
et al. (2005). The fractional uncertainties on the indi¬ 
vidual masses and radii we measure are <1%, due to the 
precise RVs that can be obtained from Keck/HIRES and 
the exquisite photometry from K2. Given that the mass 
ratio q, radius ratio r, and surface brightness ratio / are 


all nearly (but just under) unity, then the two compo¬ 
nents appear to be very nearly coeval. 

As we summarize in Table I, both stars have a to¬ 
tal broadening of their spectral lines (from rotation and 
instrumental resolution) of = 7.5 ± 0.2 km/s. The 
observations were taken in modes that produce a spec¬ 
tral resolution oi R = 36000 (FWHM) or Uy^inst ~ 3.5 
km/s, implying that the stars have an intrinsic rotation of 
<^v,rot ~ 6.6 km/s or a rotational period (given the mea¬ 
sured radii, and an assumption of spin alignment with the 
orbital plane) of Prot ~ 6 days. As can be seen in Fig¬ 
ure 3, out-of-eclipse variability seems to occur on a much 
longer timescale. This suggests that the variability might 
be a result of long-term secular changes in spot coverage 
or a beating pattern between the two stars’ complicated 
spot maps, rather than rotational modulation from a co¬ 
herent and static spot pattern. With only two intervals 
of this longer-term variation available in the K2 dataset, 
the nature of this variability remains ambiguous. 

Finally, there are now several lines of evidence demon¬ 
strating that UScoCTIO 5 is a young member of the 
Upper Scorpius OB association. As we show in Table 2 
and 4, each star has a lithium equivalent width of ^250- 
300 mA; accounting for the presence of two stars’ worth 
of continuum flux, these values are consistent with the 
500-600 mA equivalent widths presented for single M3 
stars in Upper Sco (e.g. Preibisch et al. 2001, 2002; Riz- 
zuto et al. 2015). Both stars also show Ha emission at 
levels that are consistent with the SpT-EW[Ha] sequence 
observed by Kraus et al. (2014) for non-accreting young 
stars in the Tuc-Hor moving group, as well as frequent 
flaring (Figure 3), including one small flare during a sec¬ 
ondary eclipse. 

Our measurement of the system radial velocity (7 = 
—2.64±0.07 km/s), when combined with the proper mo¬ 
tion from UCAC4 {p, = (-11.0,-18.7) ± 2.1 mas/yr; 
Zacharias et al. 2012), allows us to further test the mem¬ 
bership using kinematics. If we assume a distance of 
d = 145 ± 15 pc, then the corresponding space velocity 
with respect to the Sun is vuvw = (±2.9 ± 0.7, —14.2 ± 
2.0, —4.3± 1.4) km/s^°. The ten nearest high-confidence 
BAF members of Upper Sco have an average velocity 
of Vuvw = (±1-7 ± 0.9, -16.3 ± 0.9, -6.8 ± 0.6) km/s; 
(Rizzuto et al. 2011), while all of Upper Sco has an aver¬ 
age velocity of Vuvw = (±6.4 ± 0.5, —15.9 ± 0.7, —7.4 ± 
0.2) km/s (Chen et al. 2011); the difference between these 
measurements might point to internal kinematic sub¬ 
structure within Upper Sco. Our measurements are con¬ 
sistent with the adjacent BAF stars to within < 2 km/s 
and with the average for all of Upper Sco to within < 4 
km/s. Given the expected internal velocity dispersion of 
at least 1-2 km/s on small scales and potentially more on 
association-wide scales (e.g., Kraus & Hillenbrand 2008), 
the kinematics are therefore consistent with those ex¬ 
pected of an Upper Sco member. The proper motion 
alone is also consistent with the mean value of Upper Sco 
{p = (—9.3,—20.2) ± 0.5 mas/yr; Kraus & Hillenbrand 
2007) to within the uncertainties. We therefore further 
confirm that UScoCTIO 5 is both young and comoving 
with Upper Scorpius. 

All measurements of UVW velocities are presented in the 
sense that U is positive in the galactic anti-center direction, as 
encoded in the IDL routine gaLuvw.pro. 
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TABLE 5 

Best-Fit Model Parameters 


Model Set 

'^model 

^model 

'^best 

(7 


(Myr) 

(Mq) 

(Myr) 


BHAC15 


n 00+0.05 

D.^^_0.03 

9.7 

2.4 

SiessOO 

5.81M 

D.^U_o.o3 

10.3 

4.2 

Pisa 

6.7l?;“ 

d.^o_q Q5 

9.6 

2.4 

Padova 

IS-lta.g 

U.4i_o 05 

6.6 

1.7 

DSEP 


o O0+0.05 
D.^O_q Q4 

9.7 

2.3 


Note. — The model masses are the values for 
each star, which are identical to well within the un¬ 
certainties in the HR diagram positions. The age 
Tbest is the age along each model set’s isomass se¬ 
quence that is closest to the system’s true location 
in the HR diagram (Section 6.5); we also list the a 
level of the discrepancy. 

6. DISCUSSION 

6.1. Comparison to Stellar Evolutionary Models 

UScoCTIO 5 poses a strong challenge to the evolu¬ 
tionary models of young low-mass stars (r <1-10 Myr; 
M ~0.1-0.7 Mq), requiring predictive agreement with 
both the masses and the radii of its component stars, 
in addition to the luminosities and surface temperatures 
that are more commonly available for most young stars. 
This system also offers an intriguing older analog to 
the four previously-known eclipsing binary systems with 
young, low-mass components: 2M0535 (0.055 -I- 0.035 
Mq; Stassun et al. 2006), JW 380 (0.26 -I- 0.25 Mq; Irwin 
et al. 2007), Par 1802 (0.39-1-0.39 Mq; Gomez Maqueo 
Chew et al. 2012), and CoRoT 223992193 (0.67 -b 0.50 
Mq; Gillen et al. 2014). Those systems represent even 
more extreme youth, with inferred ages of r ~ 1-3 Myr 
for their host populations (the ONC and NGC 2264), 
but they also are complicated by the presence of circum- 
stellar disks and approach the point where models are 
distinguished more by initial conditions than by stellar 
evolutionary processes, as well as having stronger tidal 
interactions due to their shorter orbital periods. 

We specifically compare the properties of this touch¬ 
stone system (Mann et al. 2015) to the predictions of 
five sets of evolutionary models: the BHAC15 models 
(Baraffe et al. 2015), the Dartmouth Stellar Evolution 
Program or DSEP models (Dotter et al. 2008; Feiden 
et al. 2015), the Padova/PARSEC models (Bressan et al. 
2012; Chen et al. 2014), the Pisa models (Tognelli et al. 
2011), and the models of Siess et al. (2000) that while 
older, remain in common usage. In the left panels of 
Figures 9-13, we show the (nearly identical) positions of 
the two components of UScoCTIO 5 in the L-T^ff HR 
diagram, as well as the isochronal and isomass sequences 
predicted by each of the four sets of models. In each fig¬ 
ure we use the luminosity for UScoCTIO 5 calculated in 
Section 4.3 and the Teg predicted by the L03 tempera¬ 
ture scale. In Table 5, we list the corresponding model 
predictions for the mass and age of the system given its 
HR diagram position. To compute these predictions we 
adopted linear-uniform priors on both the mass (since 
the system falls near the peak of the IMF) and the age 
(implying a roughly constant star formation history in 


Upper Sco and the solar neighborhood). We also show 
the isomass sequence that the models predict for the ob¬ 
served dynamical masses in the system (M ~ 0.32Mq) 
as well as the isochronal sequence predicted by observa¬ 
tions of more massive Sco-Cen members (t ~ 11 Myr; 
Pecaut et al. 2012) that continues to gain currency over 
the longstanding predicted age determined for low mass- 
stars (t ~ 5 Myr; Preibisch et al. 2002). 

We find that none of the model isomass sequences 
successfully predict the HR diagram position of UScoC¬ 
TIO 5 to within the observational uncertainties, though 
the Padova models do predict an isomass sequence that 
agrees to within 2a. This comparison is effectively uni¬ 
variate since the isomass sequences of fully convective 
stars (which evolve along the Hayashi track) are nearly 
vertical, and hence even the Padova models still disagree 
with the observed mass at >90% in a one-sided test. In- 
triguingly, the sign and magnitude of the discrepancy is 
nearly identical between the BCAH15, DSEP, Pisa, and 
Siess models (where they under-predict the mass, as orig¬ 
inally shown by Reiners et al. 2005), whereas the Padova 
models overpredict the mass. The new Padova models 
(Chen et al. 2014) differ from the other model isochrones 
because they apply an empirical correction to the outer 
boundary correction (expressed as a Te// - t relationship 
modification from the BT-Settl photospheres) in order to 
match observed mass-radius relationships for low-mass 
dwarfs. This appears to be an over-correction in the 
case of UScoCTIO 5, as would be expected if, for exam¬ 
ple, the reason for the correction is missing atmospheric 
opacities that are most important at high gravity. How¬ 
ever, the fundamental limitations in measurement of 
(which depends on model atmosphere physics) and Ltoi 
(which depends on precise distances) will remain a limit 
on the utility of the HR diagram in testing evolutionary 
models. 

In the right panels of Figures 9-13, we show the 
(also nearly identical) positions of the two components 
of UScoCTIO 5 in the mass-radius diagram, with the 
isochronal and iso-luminosity sequences of each model 
set for comparison. These observations pose a much more 
stringent test of the evolutionary models, since both the 
mass and the radius can be determined much more pre¬ 
cisely in comparison to the dynamic range of the model 
predictions. As with the HR diagram, the two compo¬ 
nents appear very nearly coeval. The BCAH15, DSEP, 
Pisa, and Siess tracks all predict ages that are signifi¬ 
cantly younger than the newly-canonical age of t ~ 11 
Myr inferred from the upper main sequence (Pecaut et al. 
2012). As with the HR diagram, the closest agreement 
is achieved by the Padova models, which almost exactly 
reproduce the expected age. 

However, we find that none of the model sets reprodnce 
the luminosity at the give mass, with discrepancies that 
follow those of the HR diagrams. For BHAC15, DSEP, 
Pisa, and Siess, the models predict radii that are too 
small for the known luminosity and mass of the stars, in¬ 
dicating (from the Stefan-Boltzmann law) that the model 
temperatures are too high. Conversely, even though the 
Padova models predict the canonical age of r ~ 11 Myr 
for Upper Sco, the model radius for the known luminos¬ 
ity and mass is too high, indicating an underpredicted 
temperature. We therefore can demonstrate the same 
conclusions as for the HR diagrams (that there are dis- 



13 



Fig. 9. — Left: L-Tq^ HR diagram showing the measured positions of UScoCTIO 5 A+B (red error bars) and the isochronal and isomass 
sequences predicted for low-mass stars by the BHAC15 models (Baraffe et al. 2015). The components are offset slightly in Teff for clarity. 
The isomass model track corresponding to the component masses (M ~ O.32M0) and the isochrone model track corresponding to the 
currently accepted value for Upper Sco (r ~ 11 Myr; Pecaut et al. 2012) are shown in blue. Perfect agreement with the models should show 
the components of UScoCTIO 5 sitting at the intersection of the blue lines; we find that the isomass line does not match with observations, 
indicating that the Tgff predicted by the models is too high. Right: Mass-Radius diagram showing the measured positions of UScoCTIO 
5 A-l-B (red error bars) and the isochronal sequences of the BHAC15 models. As in the HR diagram, we use blue lines to show the model 
tracks for the expected isochrone (r ~ 11 Myr) and the luminosity that we measure (-L^oZ ~ O.O 66 L 0 for each star). We find that the 
position predicted by the models (at the intersection of the blue sequences) matches the mass, but not the radius; the models predict radii 
that are too small, equivalent to predicting Teff to be too high (but avoiding the systematic uncertainties of a direct comparison using Teff). 



crepancies in Teff) without relying on the systematic un¬ 
certainties of directly measuring stellar from the ob¬ 
served spectral types. 

Stars like UScoCTIO 5 remain fully convective 


throughout their evolution toward the ZAMS, and hence 
evolution on the pre-main sequence largely consists of 
dimming at constant temperature as the star contracts. 
The functional result is that luminosity is a first-order 
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Fig. 11.— As in Figure 9, but for the Padova/PARSEC models (Bressan et al. 2012; Chen et al. 2014). 



T.,i (K) 
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Fig. 12.— As in Figure 9, but for the Pisa models (Tognelli et al. 2011). 


predictor of age, while temperature is a first-order pre¬ 
dictor of mass. If a star does not fall on the appropriate 
isomass sequence for a model set, it therefore indicates ei¬ 
ther that the models are not predicting the correct value 
of Teff, or that the observational spectral types are not 
being mapped correctly to to the underlying Teff of the 
stars’ true atmospheres. If the former case is true, then it 
indicates a discrepancy in how the models handle energy 
transport. Either convection is less efficient in the inte¬ 
rior (Mullan & MacDonald 2001; Chabrier et al. 2007; 
Gennaro et al. 2012) or opacities are missing in the at¬ 


mosphere (e.g., Rajpurohit et al. 2013); either would re¬ 
sult in a smaller radius and a hotter temperature in the 
models, as we see for all but the Padova set. If the tem¬ 
perature scale is not correct (for example, due to gravity- 
dependent changes in the appearance of major molecular 
bands), then the true stellar temperature could be hotter 
or colder than predicted. 

6.2. An Empirical Constraint on the Temperature Scale 
of Young Stars 
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Fig. 13.— As in Figure 9, but for the Siess models (Siess et al. 2000). 


The model comparisons in Section 6.1 are predicated 
on a stellar effective temperature derived from the ad hoc 
young-star temperature scale of LOS. This temperature 
scale was designed to bridge the large difference between 
dwarf and giant temperatures at a given spectral type, 
and specifically to make the GG Tau quadruple and the 
low-mass stellar sequence of IC 348 appear coeval when 
compared to the models of Baraffe et al. (1998). How¬ 
ever, this temperature scale is still almost totally uncal¬ 
ibrated with observations, and there are now indications 
that the GG Tau multiple system could host additional 
components (Di Folco et al. 2014). Almost all young stars 
are too distant to measure interferometric sizes, prevent¬ 
ing a direct measurement of stellar radii. The UScoCTIO 
5 system (and other future discoveries) offer the intrigu¬ 
ing alternative of measuring geometric radii from eclips¬ 
ing binaries that are nominally model-independent. 

As we describe in Section 4.3, we estimate that the 
two components of UScoCTIO 5 have a combined lumi¬ 
nosity of Lboi ^ 0.132 ± O.O3L0. Given effective tem¬ 
peratures for each component that are virtually iden¬ 
tical (ATeff < 10 K), then the corresponding absolute 
temperatures can be derived purely from the sum of 
the radii and the Stefan-Boltzmann Law. We find that 
Teff = 0.55870 = 326OI200 "This measurement is only 
discrepant from the values predicted for dwarfs or giants 
(3100 K or 3435 K) by ±1 sigma, and hence does not 
yet provide a useful constraint on the system. However, 
the measurement is limited almost entirely by the uncer¬ 
tainty in the distance (d = 145 ± 15 pc). This system 
falls well within the brightness range where Gaia should 
deliver extremely precise parallaxes within 2 years; once 
those results are released (or if the distance can be re¬ 
fined in some other way), the system temperature can be 
described by: 


Teff = (3260lfoi^)(^)^ (1) 

6.3. An Infrared Flux Method Distance 

As we have directly measured the radii of the two 
stars in linear units, and have bolometric luminosities 
and spectroscopically estimated Tea values, we can esti¬ 
mate the angular diameters and therefore the distance 
to the system by simple trigonometry. As this is some¬ 
what dependent on reddening, we can improve this dis¬ 
tance measurement using the infrared flux method (e.g. 
Casagrande et al. 2010). We computed a model K-band 
flux for the two stars interpolating the BT-SETTL mod¬ 
els (Allard et al. 2011), using the spectroscopically de¬ 
rived effective temperatures and extinction, the K filter 
profile from Cohen et al. (2003), converting Ay to Ak 
using Savage & Mathis (1979). The final distance com¬ 
puted is 144.4±6.6pc for the L03 temperature scale or 
135.1±8.8pc for the Herczeg & Hillenbrand (2014) tem¬ 
perature scale, in good agreement with the assumed dis¬ 
tance of 145 pc to Upper Scorpius. The distance uncer¬ 
tainty is dominated by the uncertainty in the spectro¬ 
scopic effective temperature. 

6.4. The Coevality of Young Binary Systems 

The apparent coevality of UScoCTIO 5 A-l-B poses a 
counterpoint to recent suggestions that (apparent) stellar 
ages could be essentially “randomized”, which threatens 
the assumption underlying all studies of stellar evolu¬ 
tion that treat star clusters as simple stellar populations. 
Apparent non-coevality has been seen at ages of t ~ 1 
Myr for the eclipsing binary systems 2M0535 and Par 
1802 (Stassun et al. 2007; Gomez Maqueo Chew et al. 
2012), where precise characterization of the stellar pa¬ 
rameters shows that the binary components do not lie 
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on the same isochrones. Non-coevality also has been sug¬ 
gested more generally as a source of intrinsic luminosity 
spreads within binary systems and stellar populations, 
though tests of non-coevality require careful considera¬ 
tion of the non-linearity of isochrones and mass tracks, 
as noted by Gennaro et al. (2012). Two modes have been 
suggested as possible sources of this randomization, one 
tied to the assembly and one tied to the subsequent evo¬ 
lution. 

One phenomenon that could alter apparent stellar ages 
is through stochastic variations in the episodic accretion 
history of individual stars (Baraffe & Chabrier 2010), 
whereby a variable fraction of the accretion energy is 
radiated away during the accretion process and not de¬ 
posited into the star. Stars that radiate away more ac¬ 
cretion energy are left with smaller radii at a given mass 
and age. Close binaries are an imperfect test of the ac¬ 
cretion hypothesis since the stars likely accrete from a 
circumbinary disk and hence should have correlated ac¬ 
cretion histories. However, recent observations suggest 
that accretion could preferentially occur onto primary or 
secondary stars (e.g., Jensen et al. 2007) depending on 
the specific angular momentum of the accreted material, 
and hence forced-coevality might not be assured. 

The other phenomenon that could alter apparent ages 
is through variable magnetic field strengths (Mullan & 
MacDonald 2001; Chabrier et al. 2007). Strong magnetic 
fields within the stellar interior should inhibit convec¬ 
tion (preventing radial movement of charged particles), 
resulting in less efficient energy transport and a corre¬ 
spondingly larger radius and lower temperature for a 
given luminosity. Strong magnetic fields near the surface 
also should increase the starspot fraction, reducing the 
average surface temperature and hence again requiring a 
larger radius for a given luminosity. These effects could 
manifest as a correlation between fundamental properties 
(mass, radius, and temperature) and activity signatures 
such as UV or Ha emission (Stassun et al. 2012, 2014b). 

Given the very precise agreement in the apparent ages 
of UScoCTIO 5 A-l-B, combined with other internally 
coeval systems at younger ages (JW 380 in the ONC 
and CoRoT 223992193 in NGC 2264), it appears that 
ages are not significantly randomized for all young stars. 
We therefore suggest that the process forcing apparent 
non-coevality likely only occurs for a fraction of all stars. 
It has been demonstrated for visual binaries in Taurus 
(r ~ 2 Myr) that ~2/3 of all pairs appear highly coeval 
(Alog(T) < 0.16 dex or <40%; Kraus & Hillenbrand 
2009), indicating that <20% of all stars have apparent 
ages that differ substantially from the actual age. 

6.5. The Ages of UScoCTIO 5 and Upper Scorpius 

The age of the Upper Scorpius OB association, and by 
extension all young populations age-dated in a similar 
manner, is a topic of contention in the current literature. 
The ages of many populations have been set by the posi¬ 
tions of low-mass stars (such as UScoCTIO 5) in the HR 
diagram, as compared to isochrones predicted by stellar 
evolutionary models. However, the ages of upper main 
sequence stars (such as the F stars in Upper Sco) and the 
location of the age-dependent lithium depletion bound¬ 
ary seem to predict ages that are older by a factor of two. 
The most visible debate has occurred for Upper Sco it¬ 
self, with predictions of r ~ 5 Myr from low-mass stars 


(Preibisch et al. 2002; Slesnick et al. 2006) and r ~ 11 
Myr from the other methods (Pecaut et al. 2012). 

UScoCTIO 5 represents a fundamentally new data- 
point for this debate. If the evolutionary models of low- 
mass stars are indeed predicting the correct age, then 
they should also predict the correct mass and radius 
for a given luminosity and spectral type. Conversely, if 
they are predicting incorrect ages because either the lu¬ 
minosity or temperature predictions are incorrect, then 
the (time-independent) masses or the (time-dependent) 
radii might not match. This comparison is especially 
useful for low-mass stars because they evolve along the 
fully-convective Hayashi track and hence fall nearly ver¬ 
tically in the HR diagram; Teg corresponds mostly to 
mass, while L^oi corresponds mostly to age. 

As can be seen in Figures 9, 10, 12 and 13, and as 
we summarize in Table 5, the BHAC15, DSEP, Pisa, 
and Siess tracks do indeed predict a younger age for 
UScoCTIO 5 than the canonical value estimated for 
intermediate-mass stars. However, the model isomass se¬ 
quences all predict that UScoCTIO 5 should have higher 
Teff than is observed, generally approaching or exceeding 
the value that would be expected for an M4.5 giant (in¬ 
dicating that temperature scale changes alone might not 
solve this problem). The spectral line lists for low-mass 
stars are known to be incomplete, and hence the models 
must be missing opacities that would drive their pre¬ 
dicted temperatures lower. If the Teff predictions were 
shifted to bring the isomass sequences into agreement 
with the observations, given the same luminosity, then 
all three sets of models would indeed predict an age of 
r ~ 11 Myr; we list that best-fitting age along each 
model’s isomass sequence in Table 5, along with the a 
level of the discrepancy. This is unlikely to completely 
solve the problem, though, since the models can’t be 
shifted purely horizontally. Lower temperatures would 
result in less energy being radiated away, slowing the 
contraction and hence also modifying the radius and lu¬ 
minosity at a given age. 

In contrast. Figure 11 and Table 5 show that the 
Padova models predict an older age from the HR dia¬ 
gram of UScoCTIO 5, while under-predicting the tem¬ 
perature for stars that fall along its isomass sequence. 
If the model isomass sequence were shifted in Teff/Afeo/ 
space to match the observed Teff, then the inferred age 
would be younger, but still consistent with r ~ 11 Myr 
to within ^ Ict. Some authors do indeed predict that the 
young-star temperature scale should fall at even cooler 
temperatures (e.g., Herczeg & Hillenbrand 2014), which 
would bring UScoCTIO 5’s SpT-derived T^s into excel¬ 
lent agreement with the Padova grid’s existing predic¬ 
tions for the Teff and Lboi of a ~10 Myr, ~ 0.3Mq star. 

However, as the discussion above demonstrates, the 
uncertainties in UScoCTIO 5’s HR diagram position are 
dominated by the ±75 K observational uncertainty and 
the ±100 K systematic uncertainty due to the unknown 
gravity dependence of the temperature scale for young 
M-dwarf stars. A more robust comparison can be made 
in the mass-radius plane. Both quantities for the system 
can be determined with greater precision compared to 
the dynamic range of the model predictions, as well as 
without systematic uncertainty. As can be seen in the 
right-hand panels of Figures 9-13, the BHAC15, DSEP, 
Pisa, and Siess models do indeed predict ages that are 
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significantly younger than 10 Myr, though none fall as 
young as 5 Myr, while the Padova models predict the 
older age seen for intermediate-mass stars. 

The robustness of these predictions, and the direction 
of any discrepancy, can be tested by comparison to the 
iso-luminosity line for UScoCTIO 5. The intersection 
of the iso-luminosity line and the 11 Myr isochrone falls 
very nearly at the mass of UScoCTIO 5 for the BHAC15, 
DSEP, and Pisa models, at a moderately lower mass 
for the Siess tracks, and at a much higher mass for the 
Padova tracks. These results support the trend seen for 
the HR diagram - if the and radius predicted by 
the models are modified to match the isomass lines, then 
most of the tracks would indeed predict an age of r ~ 11 
Myr. We therefore conclude that while none of the model 
sets predict all system parameters perfectly, the likely 
form of the discrepancy supports the ongoing rescaling 
of pre-main sequence stellar ages in favor of older values. 

7. SUMMARY 

We have presented the discovery that UScoCTIO 5, 
a known low-mass spectroscopic binary {P = 34 days, 
M ~ O. 64 M 0 ) in the Upper Scorpius star-forming re¬ 
gion (r ~ 10 Myr), is an eclipsing system suitable 
for determination of precise stellar masses and radii. 
Based on the stellar properties {Ma = 0.329 ± O.OO 2 M 0 , 
Ra = 0.834 ± O.OO 6 R 0 , Mb = 0.317 ± O.OO 2 M 0 , Rg = 
0.810 ± O.OO 6 i? 0 ), we conclude that: 

1. There are systematic errors in the calibration 
of pre-main sequence evolutionary models. The 
BHAC15, DSEP, Pisa, and Siess models overpre¬ 
dict the Teff of young stars for a given mass by 
~200 K, or equivalently underpredict the masses 
of young stars for the given Teff by ~50%. The 
Padova models are a slightly better match, but 
are discrepant in the opposite direction, under¬ 
predicting Teff by ^100 K or overpredicting mass 
by ~25%. The discrepancies remain in the mass- 
radius-luminosity space, suggesting that the dis¬ 
crepancies likely represent intrinsic calibration is¬ 
sues rather than an uncertain temperature scale. 

2. Our geometric measurement of Tgff (derived from 
radius and luminosity) broadly agrees with the 
temperature scale for young stars, but the uncer¬ 
tainty will remain too large to refine the temper¬ 
ature scale until the luminosity can be measured 
more precisely with a Gaia distance. 


3. UScoCTIO 5 appears highly coeval, bringing the 
count among low-mass EBs to three apparently 
coeval pairs and two apparently non-coeval pairs. 
The inferred fraction of stars with spurious ages 
(~ 20 %) is consistent with the number seen for 
wider visual binaries, suggesting that processes 
which randomize apparent stellar ages do occur, 
but in a minority of cases. 


4. Taking into account the dimensions within which 
the models appear to be miscalibrated, we find that 
the age of UScoCTIO 5 appears more consistent 
with the older age of Upper Scorpius that has re¬ 
cently gained canonical status: r ~ 11 Myr. 
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APPENDIX 

The light curves of eclipsing binaries are traditionally predicted using dedicated software packages (such as the 
Wilson-Devinney code; Wilson & Devinney 1971). These packages include a large number of physical effects that can 
be relevant for some systems, such as tidal distortion, reflected light, and star spots. However, the architecture needed 
to model these effects results in a code that requires seconds to produce a single light curve, even when the features 
are turned off and even though it runs in a compiled language (Fortran). This runtime makes MGMG implementations 
onerous to execute, since long chains and multiple walkers are needed to converge and become well-mixed in the 
high-dimensional parameter space of eclipsing binary parameters. 

We therefore present here a modihcation of the formalism of Mandel & Agol (2002) that is commonly used to model 
transiting extrasolar planets. We specifically use their formalism (and their IDL implementations)^^ to predict the 
flux decrement from whichever star is being occulted, and then include an additional parameter (the surface brightness 
ratio) that, in combination with the radius ratio, can be used to add the light of the occulting star. 

http://www.astro.washington.edu/users/agol/transit.html 
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Consider two stars in an eclipsing binary, hereafter P and S, that have radii Rp and Ra (with ratio r = Rs/Rp) 
and surface brightnesses Sp and Sg at some arbitrary wavelength (with ratio s = Ss/Sp). The surface brightnesses 
will depend on TefT,p and Teff^s, but for cool stars with complicated spectra, this dependence might not be straightfor¬ 
ward. For simplicity, this algorithm fits for the ratio of surface brightnesses s, but does not attempt to calculate the 
corresponding temperatures. 

In the case where neither star is eclipsed, denoted time t = 0, the observed flux is: 

FtotiO) = FpiO) + Fa{0) = TT^Sp + TT^Sa 

Potentially also including a limb darkening term; in our case this term is identical for both stellar components, 
and hence it factors out. If the primary star is eclipsed at a time ti with an impact parameter zi(ti), then the 
Mandel-Agol algorithm can be called to calculate the fractional flux deficit with respect to the primary star’s total 
flux, fip{zi{ti)) = The other star is still contributing Fa{ti) = Fa{0), and hence the total flux observed is: 

R^ R? 

Ftot{ti) = Fp{ti) + Fa{ti) = ^J^p{zl{tl))Fp{Q) + Fs{0) = ^j,p{zi{ti))Trj^Sp + 


If we normalize this measurement by the out-of-eclipse total flux, then we can recast the observation as: 


Ftot{ti) _ + R^Sa _ fj,p{zi{ti)) + s 

FtotiO) ~ R^Sp + R^aSs “ I + r2s 


Similarly at a time t 2 during the secondary eclipse, then: 


Ri 


Ri 


Ftotit2) = Fpit2) + Fait2) = Fp{Q) + pL2{z2it2))Fai0) = '^J^Sp + /Tg (Z2 (^2 ))7r y=;|-S'' 


£12 


Ftot{t2) _ RpSp F fiait2)RaFs _ I -|-(Z2 (i2))^^S 

FtotiO) ~ Rpp + RlSa “ I + r2s 


The only input required to calculate /x(z(t)) is the impact parameter between the two stars as a function of the 
occulted star’s radius, which can be calculated analytically from the geometry of the orbit for any given orbital phase, 
and hence this formalism can be used to calculate the time-dependent fractional flux (with respect to non-eclipse 
epochs) observed for an eclipsing binary system over the course of its eclipses. Even implemented in an interpreted 
language (IDL) that is inherently much slower, this routine produces light curves 10-100 times faster than when 
using the latest Wilson-Devinney code^^ and therefore provides a suitable fast approximation when the more powerful 
features of Wilson-Devinney are not required. 
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